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The relation between the incommensurability observed in neutron scattering experiments in bilayer cuprate 
superconductors and the electronic structure is investigated. It is found that the observed incommesurability 
pattern, as well as its dependence on energy, can be well reproduced by electronic dispersions motivated by 
angle resolved photoemission data. The commensurate resonance and its contribution to the superconducting 
condensation energy are discussed in the context of these calculations. 
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I. INTRODUCTION 

One of the more controversial topics in the field of 
high temperature cuprate superconductivity is the origin 
of the incommensurability observed by neutron scatter- 
ing experiments. The original explanation of this phe.- 
nomenon was that it was due to Fermi surface nestingtl. 
Later, it was proposed that it was due to the formation 
of striped. These two explanations are so different, one 
would expect that ways of differentiating them using ex- 
isting data should be possible. In bilayer cuprates, the 
situation is even more interesting, in that a commensu- 
rate resonance is seen-ipilong with incommensurability at 
energies off resonanceHcl. 

In this paper, the Fermi surface nesting approach is 
analyzed based on tight binding energy dispersions mo- 
tivated by angle resolved photoemission spectroscopy 
(ARPES) studies, with the dynamic susceptibility cal- 
culated using the random phase approximation (RPA). 
Such tight binding dispersions are able to reproduce the 
commensurate resonance, along with the incommensura- 
bility off resonance. In particular, the incommensurabil- 
ity wavevector is found to be energy dependent, in good 
agreement with recent experimental resultsu. This would 
seem to argue against a stripes interpretation, where one 
would expect the incommensurability to be energy inde- 
pendent. Moreover, the incommensurability in the RPA 
calculations is very sensitive to the underlying electronic 
structure, which has implications for the position and 
curvature of the Fermi surface near the d-wave node. As 
for the commensurate resonance, it is a more robust fea- 
ture of these calculations, though its width in momentum 
space is sensitive to the location of tie, flat band near 
(tt, 0). In support of earlier estimatescl^LD, the change in 
exchange energy between the normal and superconduct- 
ing states in these calculations is sufficient to account for 
the superconducting condensation energy. 

In Section II, the details of the computations are given. 
In Section III, the commensurate resonance is discussed 
in relation to the electronic structure, with Section IV 
dealing with the same relation in regards to the incom- 
mensurability off resonance. Section V uses these cal- 



culations to comment on the question of the lowering of 
the exchange energy in the superconducting state. Some 
conclusions are offered in Section VI. 



II. CALCULATIONAL DETAILS 

The methodology HSed here is essentially the same as 
that of other groupsoilJ. The only difference is to use 
energy dispersions motivated by actual fits to photoemis- 
sion measurements. As in these earlier studies, one first 
determines the non- interacting susceptibility, which., i 
the superconducting state is composed of three termstd 
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where efc is the dispersion, Aj, the superconducting gap, 
Ek the quasiparticle energies {^/e^'+'K^), and / the 
Fermi function. The three terms are due to quasiparti- 
cle scattering, quasiparticle pair annihilation, and quasi- 
particle pair creation, respectively. At low tempera- 
tures, only the last process contributes. Note this equa- 
tion makes a severe approximation of treating the single 
particle spectral function as a delta function. On the 
other hand, at low temperatures, photoemission data in- 
dicate ouasiparticle peaks at all k vectors near the Fermi 
surfacell3 (though alonaihe d-wave node direction, there 
is still some controversyllJ) . As the incoherent part of the 
spectral function is unlikely to lead to sharp structure in 
q and these earlier studies based on Eq. 1 are followed, 
but the reader should keep in mind that there are cer- 
tainly quantitative, and perhaps qualitative, differences 
between this "quasiparticle" susceptibility and the true 
one which warrant future investigation. 

To evalulate Eq. 1, i5 is replaced by iT in the energy 
denominators. The resulting well behaved function is 
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then summed over a 400 by 400 mesh in the Brillouin 
zone, r, which is a crude representation of broadening 
due to interactions and disorder, was taken to be 2.4 
meV. SmaUer values of F simply lead to sharper structure 
in xq which requires a denser k mesh in the sum. 
The interacting susceptibility in the RPA is given by 



1 - J{q)xo{q,^) 



(2) 



were J{q) is the spin-spin response function. Many sLud,- 
ies set J{q) to a constant, J, but in several treatmentsSO, 
J{q) is assumed to be of the form Jq = — J(cos(qj,a) + 
cos(gj,a))/2 due to superexchange between near-neighbor 
Cu sites, so calculations were performed with both 
J{q) = J and J{q) = Jq. The use of Jq tends to suppress 
incommensurability since Jq is largest at the commensu- 
rate wavevector (7r,7r). 

For energy dispersions; a uumber of model dispersions 
present in the literaturcE3E30'EII were analyzed. These 
are based on tight binding functions, and in the more 
sophisticated models designed to fit photoemission data, 
they contai]p_aip_fixpansion up to a real space lattice vec^ 
tor of (2,2)y"El. The first to be considered is oneN 
which does a very good job of modeling the norrpil state 
ARPES dispersion in BiaS'raCaCuaOs (612212)113. This 
is reproduced in Table F Note that this model dispersion 
has no bilayer splitting, consistept with its lack of ex- 
perimental observation in Bi2212lia. J is then chosen to 
yield a resonance at a particular energy. A d-wave energy 
gap is assumed of the form = A{cos{kx) — cos{ky))/2. 
An s-wave gap strongly reduces the tendency to obtain a 
commensurate pesonance because of the BCS coherence 
factors in Eq. llij'E3. 



III. COMMENSURATE RESONANCE 

The first question addressed concerns the com- 
mensurajtej j cson ance, first observed in Y Ba2Cu30i 
(YBCO)By^By, and recently in Bi2212yy. The con- 
ditions under which a commensurate response is obtained 
for the non-interacting susceptibility in a simple t, t' tight 
binding model is given in Rcf. 10, and are roughly obeyed 
for the more sophisticated dispersions considered here. 
For simplicity, consider first the normal state. Then, the 
condition reduces to whether the Fermi surface is cen- 
tered at (0,0) or at (tt, tt). In the former case, the re- 
sponse becomes incommensurate, basically because the 
Fermi surface and its shadow image, displaced by a 
wavevector (tt, tt), no longer intersect (the intersection 
points are refered to as hot spots in the literature) . Once 
a superconducting gap opens, jtiiere is a correction to this 
condition which scales like A^t2l. For the dispersion ana- 
lyzed here, this correction corresponds to about 5 meV. 
That is, once the band at (tt, 0) becomes more than 5 
meV above the Fermi energy, the response becomes com- 
pletely incommensurate. All models for YBCO, and most 



models for Bi2212, are in favor of Fermi surfaces centered 
about (tTjTt). Aft-exception is a recent model by Chuang 
et al for Bi2212ca. In this model, the Fermi crossing is 
far enough away from (tt, O) along the (0, 0) — (tt, 0) line 
that an incommesuratc response is expected, unlike what 
is observed experimentallycj. 

The non-interacting susceptibility having its maximum 
response at (tt, tt) leads to a maximum at (tt, tt) in the 
interacting susceptibility as well. For a d-wave gap, the 
coherence factors in Eq. 1 are maximal on the Fermi sur- 
face, due the sign change of the gap under (tt, tt) trans- 
lation. This leads to an abrupt rise in Imxo above some 
threshold energy value (twice the superconducting energy 
gap at the hot spots), which is not present for the s-wave 
case. By Kramers-Krottig transformation, the abrupt rise 
causes a peak in Rexo^- Therefore, if J is large enough, 
a zero will occur in the real part of the denominator of 
Eq. 2 for energies smaller than the threshold. As Im^o is 
small below the threshold, this leads to a resonance like 
behavior in Im^. That is, in general, one expects a reso- 
nance at a energy smaller than 2lS.max- These points are 
illustrated in Fig. 1. Note that all the peak structures in 
this plot become more pronounced if the broadening, F, 
is reduced in size. 

An interesting evolution of this resonance occurs as the 
chemical potential, /i, is varied. Using the same model 
dispersion, J was adjusted for each /i in such a way that 
the maximum in Imx occurs at 39 meV (with Amaa: taken 
to be 35 meV). In Fig. 2, a 2D contour plot is shown at 
the resonance energy for various positions of the band 
at (7r,0), e(7r,0). As e(7r, 0) approaches from below, 
the resonance broadens in momentum space, and even- 
tually transforms from a circular pattern to a square-like 
pattern. As the band crosses /i, the resonance pattern ro- 
tates to a diamond shape, and then far enough beyond, 
a completely incommensurate response is achieved. In 
agreement, with the analytic results on the t,t' model 
for Imxo^, the value of e(7r, 0) relative to /i where the 
crossover to incommensurate behavior occurs is found to 
be proportional to A^. The width of the resonance in 
momentum space is also sensitive to other factors. The 
above calculations assumed J{q) — J. If J{q) — Jq is 
used instead, the resonance would become narrower in 
q, which is obvious from the q dependence of Jq, Jq be- 
ing maximal at (7r,7r). Also, the dispersion used in the 
above calculation is based on normal state ARPES data. 
In the superconducting state, though, it is, known that 
the dispersion becomes fiatter near (tt, 0)E3, which acts 
to broaden the resonance in momentum space. From ex- 
periment, the resonance inj^i2212 is broader in momen- 
tum space than in YBCOEj. In the context here, this 
would imply that for Bi2212 relative to YBCO: (1) the 
band at (tt, 0) is closer to the chemical potential, (2) the 
dispersion at (tt, 0) is flatter, and/or (3) J{q) is fiatter 
near (7r,7r). A quantitative comparison of experiment to 
theory will be discussed in Section V. 

At this stage, nothing has been said about the qz de- 
pendence of the resonance. Experimentally, the reso- 
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nance only appears in the odd {qz = tt) channel. There 
are three ways this could occur in the context of the RPA 
calculations. First, J{q) could be larger in the odd chan- 
nel than the even one due to interplane exchangaj. In this 
case, the pronounced gap in the even channel would be 
associated with the threshold energy discussed earlier for 
Imxo (twice the energy gap at the hot spots). Altbpugh 
experimentally, J{q) is larger in the odd channelEj, the 
value of the interplane exchange integral (Jj_) is small 
enough that there would be qualitative problems with fit- 
ting the even channel data. For instance, for both disper- 
sions listed in Table I, given a value of J needed to obtain 
a resonance at 39 meV in the odd channel, a 20% reduc- 
tion in J to simulate the even channel still results in a 
pronounced resonance. An exception is the dispersion of 
Ref. where the same analysis leads to no resonance in 
the even channel. This occurs because the peak in Rexo 
in Fig. 1 is very shallow. For the same reason, though, a 
value of J much larger than experiment is needed (~ 1 
eV) to obtain a resonance in the odd channel. 

The second way would be to recognize that the odd 
channel corresponds to connecting bonding to antibond- 
ing states in Eq. 1, whereas the even channel conuects 
bonding to bonding and antibonding to antibondingB. In 
this context, bonding and antibonding refer to the bilayer 
splitting of the electronic structure. As mentioned above, 
there is little evidence for such splitting from AEPES 
data in Bi2212, even in the superconducting statella (al- 
though recently, this result has been challengedEj) . Still, 
several calculations including bilayer splitting were an- 
alyzed. In the first case, dispersion two in Table I 
was taken to be the bonding band, and the antibond- 
ing dispersion was gotten by adding a constant in such 
a way that its Fermi surface crossing occured along the 
(0, 0) — (tt, 0) line. In the second case, the bilayer split dis- 
persion of Ref. |l^ was used. In the first case, there was 
virtually no change in the results. In the second case, 
there were some minor quantitative differences (i.e., a 
pronounced resonance still existed in the even channel). 
It should be remarked that although the evidencej_£pr bi- 
layer splitting is considerably greater in YBCOEH, the 
interpretation of ARPES data in this case is more con- 
troversial due to the contribution from the chains, as well 
as surface related problems. 

The third possiblity is that the (AfcAfe+g) correlator 
in Eq. 1 only contributes to the odd channel for some 
reason not apparent at the moipentEil. In this regard, it 
is interesting to note that JankoEd has recently predicted, 
based on thermodynamic arguments, that the resonance 
will be strongly suppressed for fields along the c-axis. As 
Janko discusses, the likely source of this field dependence 
is that the (AkAk+g) correlator in Eq. 1 is sensitive to 
phase coherence. This has also been suggested to be the 
case in the context ofjnterpreting c-axis optical conduc- 
tivity measurements^ This same phase coherence sensi- 
tivity might be linked to why the resonance only appears 
in the odd channel. 



IV. INCOMMENSURATE RESPONSE 

The next question concerns the incommensurate be- 
havior. jOfF resonance. Only in YBCO is this known in 
detaiflp. Incommensurate behavior is observed in Bi2212 
as wellcil, but the 2D pattern of this in q space, as well 
as its energy dependence, is unknown at this time. Inter- 
estingly, in YBCp, the incommensurate wavevector de- 
pends on energyQ. It basically has an "hourglass" shape, 
with the incommensurability wavevector approaching the 
commensurate value as the resonance is approached from 
below, then again splitting out above the resonance en- 
ergy. In the RPA calculations, although incommensura- 
bility above the resonance energy is a robust result, being 
related to the commensurate/incommensurate discussion 
of the previous Section, the incommensurability below 
the resonance energy is a different matter altogether. In- 
commensurability belew resonance has been present in 
previous calculationaSJ, and the explanation for it within 
the RPA contest was given in a recent paper by Brinck- 
mann and LealU, with the calculations presented here 
in agreement with their picture. In the superconduct- 
ing state, the constant energy contours at low energies 
are very elongated due to the large ratio of the Fermi 
velocity at the d-wave node to the slope of the supercon- 
ducting energy gap at the node. This velocity ratio has 
been exaerimentally determined by ARPES to be 20 for 
Bi2212Ea, and the same value has been extracted froiii 
low temperature thermal conductivity measurementsLJ. 
The latter measurements have also determined the ve- 
locity ratip to be 14 in YBCO. As shown by Brinckmann 
and Leetll, the incommensurability is due to nesting be- 
tween the energy contour about the d-wave node to the 
same contour displaced by a wavevector Q = (tt, tt), that 
is, Ek ~ w/2, Ek+Q '~ '^/2, where lo is the neutron en- 
ergy. The wavevector is incommensurate since the Fermi 
surface at the node is displaced away from the {n/2, tt/2) 
points. 

For the dispersion used in Fig. 2, though, no in- 
commensurability is found below the resonance energy, 
though incommensurability does occur above the reso- 
nance energy. This can be traced to the fact that for 
this dispersion, the Fermi surface is too curved at the 
node. This leads to low energy contours which have a 
"banana" shape, thus destroying nesting, as illustrated 
in Fig. 3. Moreover, for this dispersion, the Fermi cross- 
ing at the d-wave node corresponds to a wavevector of 
(0.37, 0.37)7r, which would yield a larger incommensura- 
bility than is typically observed. Both of these problems 
can be corrected if the Fermi wavevector is pushed closer 
to the (7r/2,7r/2) point. This has the effect of reducing 
the Fermi surface curvature, thus enhancing nesting, and 
also reducing the magnitude of the incommensurability. 
Reducing the Fermi velocity also aids the nesting, but 
this is a smaller effect. 

To investigate these points further, several modifica- 
tions to the dispersion used in Fig. 2ll3 were made. First, 
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the very flat dispersion of the superconducting quasipar- 
ticle states near (tt, 0) was incorporated by setting e(7r, 0) 
to -5 meV, and then invoking the condition that the cur- 
vature of Efe is zero along the and ky axes at this 
point. To obtain an incommensurability more relevant 
to experimental data, the Fermi wavevector at the node 
was pushed out to (0.41, 0.41)7r. At the same time, the 
Fermi velocity was reduced from 1.6 eVA to 1.0 eVA, 
but as discussed above, this has a smaller influence on 
the results. This dispersion (two) is also listed in Table 
I, with the resulting Fermi surface and low energy con- 
tours shown in Fig. 3 as well. Note the flatter energy 
contours as compared to the previous dispersion. 

In Fig. 4 (left panels), 2D contour plots are shown for 
this dispersion at three different energies: on resonance, 
and 10 meV above and below resonance. Note that the 
incommensurability pattern below resonance has a strik- 
ing "baseball diamond" shape, very similar to what has 
been observed experimentally in YBCOH. This pattern 
was generated assuming J{q) — J. If J{q) — Jq is used 
instead (right panels) , another maximum develops at the 
commensurate wavevector, and is or is not the global 
maximum depending on the particular energy. In this 
context, a number of model dispersions were analyzed. 
Some exhibit completely commensurate behavior below 
resonance, others completely incommensurate behavior, 
and others still mixed behavior. That is, in the RPA 
context, the incommensurability below resonance is very 
sensitive to the electronic strucuture, as well as to the 
momentum dependence of J{q). 

In Fig. 5 (left panel), the wavevector along the (tt, tt) 
direction where Imx is maximal is plotted versus energy, 
using the same parameters as Fig. 4. Note the distinct 
"hourglass" shape of the pattern, which has recently been 
observed experimentally in YBCOa. The striking agree- 
ment of this pattern with experiment is a strong argu- 
ment in favor of an RPA-like interpretation of the data. 
On the other hand, the RPA calculations do suffer from 
some quantitative problems. In the right panel of Fig. 5, 
the intensity at the wavevector where Imx is maximal 
is plotted versus energy. Note the extremely sharp drop 
as one moves ofL resonance. Experimentally, this drop is 
less pronouncedEl. 

As suggested above, in the RPA context, the neutron 
scattering results are a sensitive probe of the electronic 
structure. This raises the question of whether the in- 
commensurability structure is in quantitative agreement 
with ARPES results or not. Unfortunately, all neutron 
scattering studies on this issue but one have been done 
on YBCO. ARPES results on YBCO are stiU somewhat 
controversial because of surface related issues not present 
in Bi2212. For instance, the tight binding dispersions 
proposed in Ref. |l^ do not support incommensurate be- 
havior below the resonance energy, again because of too 
strong a curvature of the Fermi surface around the node 
(that is, the low energy contours are too curved to sup- 
port nesting). Moreover, published ARPES results on 
YBCO indicate a Fermi crossing along the node which 



would result in an incommensurate wavevector which is 
too far displaced from (7r,7r) relative to experiment. It 
would, of course, be desirable in YBCO to exploit the ad- 
vent of the Scienta high momentum resolution detectors 
to revisit this issue. 

In Bi2212, the ARPES Fermi surface (in the vicinity 
of the d-wave node, at least) is better agreed upon. As 
for neutron results, only one study has been offeredtj. 
In that experiment, a rod of crystallites aligned along 
the (1, 1) direction was measured. This indicated an in- 
commensurate wavevector at low temperatures of around 
(0.82, 0.82)7r. Note that the two dimensional pattern 
of the structure is not known from these data. More- 
over, energy information is also not known, so there is 
always the possibility that the incommensurability being 
observed is above resonance. Lacking further data at this 
stage, let us assume that incommensurability is indeed 
being observed below resonance. Then, for the ARPES 
data to be consistent with the neutron wavevector, the 
Fermi crossing along the node would have to be at about 
(0.40, 0.40)7r. This is close to a recentlj: reported value of 
(0.39, 0.39)7r using a Scienta detectortil. So, the reported 
wavevectors from the ARPES and neutron experiments 
are certainly within current error bars, and therefore con- 
sistent at this stage. It would be highly desirable to: (1) 
have neutron data on the incommensurability on single 
crystal samples of Bi2212, and (2) to have more accu- 
rate measurements of both the Fermi wavevector, and 
the curvature of the Fermi surface around the node, from 
ARPES data, to further explore this point. 

V. SUM RULE AND CONDENSATION ENERGY 

Finally, the intriguing question of the relation of the 
commensurate resonance to the superconducting conden- 
sation energy can be treated very straightforwardly in 
the RPA context. The advantage of these calculations is 
that all wavevectors and energies are accounted for, and 
therefore these calculations provide an important check 
to the ideas proposed in Refs. H and ^. 

To begin with, x has to be converted to units quoted 
in neutron scatterering work. This is achieved by mul- 
tiplying Eq. 2 by the appropriate matrix element. For 
simplicity, consider the zz matrix element, which is 
E.g^t^WS.Wr. For 5 = 2, 5 = 1/2, this reduce^ 
to 2/1^, as do the other two {xx, yy). The sum rulecZI 
can now be checked by summing Eq. 2 over the zone, 
integrating uj out to the band edges, and multiplying by 
2//^. The dispersion used in Fig. 5 is employed here, 
with a resonance at 39 meV obtained by setting J to 159 
meV, a value comparable to experimental estimates of J 
in YBCOEj. The normal state is calculated by simply 
setting A to zero. This is a gross approximation, since it 
assumes quasiparticle states at T=0 for the normal state, 
which is highly improbable. In principle, a more accurate 
represenation of the normal state could be simulated by 
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increasing F, but for simplicity, this is not done here. For 
J{q) = Jq, the sum is in the normal state, 1.64/i^ 

in the superconducting state. For J{q) = J, the sum is 
in the normal state, 1.87/i^ in the superconduct- 
ing state. So, to within a few percent, the sum rule is 
satisfied by the RPA calculations. For a local mo men t 
system, we would expect the value irg'^ figS{S 
which for g = 2, S — 1/2 reduces to n. That is, the above 
values range from 52% to 61% of the local moment result. 
This reduction is to be expected, since Eq. 2 is based on 
itinerant electrons. 

A useful comparison to experiment is to integrate 
Imx(7r,7r) over energy. Restricting to a 50 meV energy 
range, a value of 1.9^% per plane is calculated, com- 
pared to an expepmental value of 0.95//^ in Bi2212 aad, 
0.8^^ in YBCOfea. This somewhat large overshootES 
is reduced when looking at the local susceptibility (q 
integrated). The maximum (at the resonance energy) 
per plane from the above calculation is 9.8/i^/eV for 
J{q) = Jq, 14.9/^^/eV for J{q) — J. This is comparable 
to the 9.5^^/eV value quoted for underdoped YBCOB. 
10 meV below resonance, the numbers are 1.1 and 2.2, 
respectively, compared to an experimental value of 2.50. 
As noted in the previous Section, for any given calcula- 
tion, the intensity appears to drop off below resonance 
faster than experiment. Another useful comparison is to 
look at the full width half maximum of the resonance in 
q space. The calculated values are 0.34A^^ for J{q) = J, 
0.23A~^ for J{q) = Jq. This is to be compared to ex- 
perimci^al values of 0.52l-i for Bi2212 and 0.25A~^ for 

YBcoy. 

As noted in earlier workETLl, the exchange energy con- 
tribution to the free energy in the t-J model, denoted 
as Ex, is obtained by multipliying Imx as defined in 
Eqs. 1 and 2 by — 3Jq/27r, integrating over energy, and 
summing over the zone. Its contribution to the conden- 
sation energy is E^^ — E^ , where NS represents an 
extrapolation of the normal state to zero temperature. 
For the case considered above, a value per plane of 28K 
is found if J{q) = Jq is used, 59K if J[q) = J. That 
is, within the t-J model context, the exchange energy is 
indeed lowered in the superconducting state, with the 
calculated-values somewhat larger than those based on 
the dataoEI. The advantage of the current calculation is 
that this difference can be looked at as a function of q 
and uj. In Fig. 6a, the zone sum of Jqlmxq is plotted as 
a function of w [J{q) = J) for normal and superconduct- 
ing states. Note that the two merge at 2lS.rnax: that is, 
the energy difference is confined to energies below this. 
In Fig. 6b, the u integrated quantity is plotted versus q. 
Note the contribution near (tt, tt) due to the resonance. 
In addition, a normal state contribution at low q is re- 
moved in the superconducting state. If J{q) = Jq is used 
instead (Fig. 6d), the low q structure in the normal state 
is suppressed, which is why the condensation energy is 
half that of the J{q) — J case. This additional low q 
structure is likely not relevant to experiment since there 



is no evidence for it, and it is doubtful whether the true 
J{q) is large at small wavevectors. 

The implications of Fig. 6 is clear. The dominant con- 
tribution to the exchange energy part of the condensa- 
tion energy is due, to the resonance. This in aiupport 
of previous workHfl. Moreover, as noted beforeOlII, the 
estimated value, such as from the above calculation, is 
more than sufhcient to account for experiment, as the 
total condensation energy has been estimated to be only 
3K per copper oxide piane from specific heat data for 
optimal doped YBCOEj. (It should be mentioned that 
within the t-J context, if the exchange energy is lowered 
in the superconducting state, thsH it is expected that 
the kinetic energy would increaseeHl. ) It is also interest- 
ing to remark that the RPA calculations appear to be in 
greater agreement with the idea of Ref. ^ (that the res- 
onance dominates the exchange energy difference) than 
related calculations based on the spin-fermion modelcil, 
despite qualitatively similar physics. In the latter model, 
the merger noted in connection with Fig. 6a occurs at 
much higher energies (of order J). 

Finally, some comments in regards to the nature of the 
resonance mode are in order. The quantum numbers of 
the resonance correspond to an excited triplet {S = 1) 
pair with center-of-mass momentum Q = {tt,tt), since 
the BCS ground state is S" = 0, Q = 0, and the reso- 
nance is seen by spin flip scattering. On the other hand, 
the question of whether the mode is actually a particle- 
particle mode is a more delicate question, due to particle- 
hole mixing in the superconducting stateEj. The RPA 
calculation assumes that the underlying action is in the 
particle-hole channel. Is there any experimental support 
for this? As Demler and Zhang point outc3, this question 
can only be indirectly answered by neutron scattering, 
as the neutrons only couple to the particle-hole channel. 
Their argument is that since the resonance only appears 
below Tc, and since particle- hole mixing occurs below Tc, 
then a particle-particle mode would only become visible 
in neutron scattering below Tc, in agreement with experi- 
ment. Of course, in the RPA calculations presented here, 
the mode also appears only in the superconducting state, 
since it is a consequence of the BCS coherence factors in 
Eq. 1 (that is, the mode is best not thought of as just a 
spin- wave mode). Now, although it is true that a spec- 
tral gap opens up above Tc in underdoped materials (the 
pseudogap) , the RPA calculations are basied on quasipar- 
ticle states, which only appear below TJ^. So, in that 
sense, the Demler-Zhang argument does not necessarily 
resolve this issue. 

On the other hand, angle resolved photoemission does 
not have the same restriction as neutron scattering. 
Strong arguments have been made that the the dramatic 
change in the spectral lineshape observed in such data be- 
low Tc are a direct consequence of thfi inte raction of the 
electrons with the resonance modeES'Cj'EZl. This comes 
from the dressing of the electron propagator by the res- 
onance. Note that unlike neutrons, there is nothing re- 
stricting this coupling to be in the particle-hole chan- 



5 



nel. That is, one expects that if the mode were particle- 
particle in nature, the dominant coupling to the self- 
energy of the electrons would be in the particle-particle 
channel. Let us now think about the simple limit that 
the mode energy goes to zero. Then, the dispersion of the 
higher binding energy feature seen in ARPES data (the 
so-called hump) will have a dispersion given by solving a 
very simple 2 by 2 secular equatiorO. In the particle-hole 
case, the diagonal elements will be Cfc and e/c+g, where 
Q = (tt, tt). In the particle-particle case, the second el- 
ement would become —€-k+Q instead. To first approxi- 
mation in both cases, the off-diagonal elements, denoted 
as At/, are taken to be constants. In Fig. 7, the dis- 
persions obtained from both secular equations are plot- 
ted using dispersion one of Table I for e/c, for a typical 
value of /S.U (100 meV). For the particle-hole case, the 
resulting dispersiojijs very similar to what is observed 
in photoemissioneSa. But, the particle-particle case has 
no resemblance at all to the data. 

Now, this argument does not definitively rule out a 
particle-particle explanation for the mode, since it is con- 
ceivable that the dominant coupling of the mode to the 
electrons could still be in the particle-hole channel be- 
cause of the interaction vertices. Still, the above argu- 
ment is certainly very suggestive of a particle-hole origin 
for the mode, and would also more naturally explain how 
the higher binding energy feature crosses ovfj)' t|0 the Mott 
insulating gap as the doping is reducedr^^' 



VI. CONCLUSIONS 

In summary, this paper has shown that the RPA treat- 
ment of the dynamic susceptibility gives very useful in- 
sight into neutron scattering data in the cuprate super- 
conductors. In such a framework, the neutron data are 
a sensitive probe of the underlying electronic structure. 
Using dispersions motivated by angle resolved photoemis- 
sion data, a natural explanation is found for the magnetic 
resonance observed by neutrons, as well as the incom- 
mensurability seen off resonance. Moreover, ttwse calcu- 
lations are in support of previous suggestiona3i3 that the 
resonance mode provides the dominant contribution to 
the change in the exchange energy between the normal 
and superconducting states. The advantage of the RPA 
calculations is that they provide quantitative information 
on all of these issues. In particular, the current study sug- 
gests that the two momentum resolved probes used for 
the cuprates, ARPES and neutrons, are strongly related 
to one another, and are consistent within current experi- 
mental error bars. With the advent of higher momentum 
resolution detectors in ARPES, and large enough sam- 
ples for neutron studies, this connection in the future 
can be studied with much greater precision, especially in 
the case of Bi2212. 
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TABLE I. Tight binding dispersions based on fitting 
ARPES data. The first two columns list the coefficient, Ci, of 
each term (eV), that is e(fc) = ^ CiTpXk), with "one" a previ- 
ous fit to normal state ARPES datatJ, and "two" a modified 
fit as discussed in the text. The last column lists the basis 
functions (the lattice constant a is set to unity). 



one 



two 



0.1305 
-0.5951 

0.1636 
-0.0519 
-0.1117 

0.0510 



0.0879 
-0.5547 

0.1327 

0.0132 
-0.1849 

0.0265 



i 

^(cos kx + cos ky) 
cos kx cos ky 
I (cos 2kx + cos 2ky) 
i (cos 2kx cos ky + cos kx cos 2ky 
cos 2kx cos 2ky 
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FIG. 1. Imxo, Rexo, and Imx at g = (tt, tt) using dispersion 
one listed in Table I, with J=0.3 eV, r=2.4meV, A=35meV, 
and r=13K. Imx units (Eq. 1) are states/eV/CuO plane, 
which applies for all figures. 
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FIG. 2. Dependence of Imx(5x,9!/) at 01=39 meV on 
e{n, 0) - n: (a) -34 meV, (b) -10 meV, (c) meV, (d) +5 

mcV, (c) +10 mcV, and (f) +20 moV (with (a) the value 
from dispersion one of Tabic I). Same parameters as Fig. 1, 
except that for each plot, J{q) = J has been adjusted so that 
the maximum in Imx is at this u). Note that qx and qy are in 
units of TT. 
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FIG. 3. Fermi surface (dashed curve) and superconduct- 
ing state energy contours {Ek,Ek+Q = ijj/2, u}=29 meV), for 
(a) dispersion one and (b) dispersion two in Table I. Note 
Q = (1, 1) in these units. 



1.50 



0.50 



1 .50 




I ' I ' I r 

0.50 0.75 1.00 1.25 1.50 
qx 



0.50 



I ' I ' I ' 
0.50 0.75 1.00 1.25 1.50 
qx 




T ' \ ' r 

0.50 0.75 1.00 1.25 1.50 
qx 



1 ' I ' r 

0.50 0.75 1.00 1.25 1.50 
qx 



FIG. 4. Imx{qx,<ly) at three energies: (a) 29 meV, (b) 39 
meV (resonance) , and (c) 49 meV, for dispersion two of Table 
I, with J=159 meV {J{q) = J), (d), (e), and (f) are the same, 
but with J{q) = Jq. 
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FIG. 5. (a) q vector along (tt, tt) direction where Imx(g, q) 
is maximal versus lu. (b) The magnitude of Imx in (a) versus 
uj. Same parameters as Fig. 4, with J{q) = J. 
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FIG. 6. (a) Jqlmx{q,uj) and (b) J^^"""^"^^ JJmx(g,t^) 
for J(q) = J. Same parameters as in Fig. 5. (c) and (d) are 
the same, but with J{q) — Jq. NS is the normal state, SC 
the superconducting state. The exchange energy contribution 
to the condensation energy would be obtained by integrating 
(a) or (c) over uj, multiplying by — 3/(27r), and subtracting 
SC from NS. Note (b) and (d) are in eV units. 
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(0,0) (7t,7t) (7t,0) (0,0) 

FIG. 7. Solution of particle-hole (ph) secular equation, 
with diagonal elements tk^tk+q^ and particle-particle (pp) 
secular equation, with diagonal elements tk,—e-k+Q, with 
Q = (tt, tt). In both cases, the off-diagonal elements were 
taken to be A(7=100 meV. Dispersion one of Table I was em- 
ployed (ai=0). 
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